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The dark energy universe equation of state (EOS) with inhomogeneous, Hubble parameter de- 
pendent term is considered. The motivation to introduce such a term comes from time-dependent 
viscosity considerations and modifications of general relativity. For several explicit examples of such 
EOS it is demonstrated how the type of future singularity changes, how the phantom epoch emerges 
and how crossing of phantom barrier occurs. Similar cosmological regimes are considered for the 
universe with two interacting fluids and for universe with implicit EOS. For instance, the crossing 
of phantom barrier is realized in easier way, thanks to the presence of inhomogeneous term. The 
thermodynamical dark energy model is presented where the universe entropy may be positive even 
at phantom era as a result of crossing of w = — 1 barrier. 

PACS numbers: 98.70.Vc 
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I. INTRODUCTION 

The increasing number of evidences from the observa- 
tional data indicates that current universe lives in a nar- 
row strip near w = — 1 (where w is the equation of state 
(EOS) parameter), quite probably being below —1 in so- 
called phantom region. It is also assumed that modern 
universe is filled with some mysterious, negative pressure 
fluid (dark energy) which represents about 70 percents of 
total energy in the universe. (The simplest, phcnomcno- 
logical approach is to consider that this fluid satisfies 
to EOS with constant w). The origin of this dark en- 
ergy is really dark: the proposed explanations vary from 
the modifications of gravity to the introduction of new 
fields (scalars, spinors, etc) with really strange proper- 
ties. Moreover, forgetting for the moment about the ori- 
gin of dark energy, even more-less satisfactory mechanism 
of evolving dark energy is missing, so far. At best, each 
of existing theoretical models for dark energy explains 
some specific element(s) of late-time evolution, lacking 
the complete understanding. Definitely, the situation 
may be improved with the new generation of observa- 
tional data when they will present the realistic evolving 
EOS of dark energy during sufficiently large period. 

The most strange (if realistic) era in the universe evo- 
lution is phantom era. There are many attempts to de- 
scribe the phantom cosmology (see, for instance, Q, Q 
and references therein), especially near to future, finite- 
time singularity (Big Rip) which is the essential element 
of classical phantom cosmology. (Note that quantum ef- 
fects may basically provide the escape from future, finite 
type singularity, for recent discussion, see0,Q). Unfor- 
tunately, the easiest way to describe the phantom cos- 
mology in the Lagrangian formulation leads to the ne- 
cessity of the introduction of not very appreciated scalar 
with negative kinetic energy |5j- Another, easy way is 
to use some phenomenological EOS which may produce 
dark epoch of the universe (whatever it is). It is re- 



markable that such description shows the possibility of 
other types of future, finite type singularity. For in- 
stance, even when EOS is suddenly phantomic (near to 
rip time where negative pressure diverges), the sudden 
singularity occurs There may exist future singulari- 
ties where energy/pressure is finite at rip time, for clas- 
sification of future singularities, see They may occur 
even in modified gravity at late times, see[7| for explicit 
examples. Nevertheless, it is remarkable that effective 
phantom phase may be produced also in string-inspired 
gravities |g. 

The present paper is devoted to study the phantom 
cosmology and related regimes (for instance, crossing 
of phantom divide) when phenomenological equation of 
state of the universe is inhomogeneous. In other words, it 
contains terms dependent explicitly from Hubble param- 
eter (or, even from its derivatives). Definitely, one needs 
quite strong motivation for such modification of dark en- 
ergy EOS. The first one comes from the consideration 
of time-dependent bulk viscosity 0, 0- (For earlier dis- 
cussion of cosmology with time-dependent bulk viscosity, 
see see also 11].) Actually, it was constructed the specific 
model of dark energy with possibility of crossing of phan- 
tom divide due to time-dependent bulk viscosity . The 
construction of EOS from symmetry considerations 0] 
indicates to the necessity of some inhomogeneous correc- 
tion. Finally, big number of gravities: from low-energy 
string effective actions to gravity with higher derivative 
terms or with inverse terms on curvatures modifies the 
FRW equations in requested form. 

The paper is organized as follows. In the next sec- 
tion we consider spatially-flat FRW universe filled by the 
ideal fluid with specific, dark energy EOS 0. Short 
review of four types of future singularity for different 
choices of EOS parameters is given, following to ref.Q. 
The inhomogeneous term of specific form is introduced 
to EOS. The role of such term in the transition of differ- 
ent types of singularity to another ones is investigated. 
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The cosmological regimes crossing phantom barrier due 
to such terms are explicitly constructed. Finally, the de- 
pendence of the inhomogeneous term from Hubble pa- 
rameter derivatives is briefly discussed as well as emerg- 
ing oscillating universe. Section three is devoted to the 
study of similar questions when FRW universe is filled 
by the interacting mixture of two fluids. The modifica- 
tion of two fluids EOS by inhomogeneous term is again 
considered. The explicit example of late-time cosmology 
(which may be oscillating one) quite naturally crossing 
the phantom divide in such a universe is presented. It 
is interesting that inhomogeneous term may effectively 
compensate the interaction between two fluids. In the 
section four we discuss the FRW cosmology admitting 
the crossing of barrier w = — 1 due to specific form of the 
implicit dark energy EOS proposed in ref.|l3(. Again, the 
generalized, Hubble parameter dependent EOS is consid- 
ered. Some thcrmodynamical dark energy model passing 
the barrier w = — 1 is constructed, based on above EOS. 
It is demonstrated that in such a model the universe 
entropy may be positive even during the phantom era. 
Some summary and outlook are given in the discussion 
section. The Appendix deals with couple simple versions 
of modified gravity which may predict the requested gen- 
eralization of EOS. 



II. FRW COSMOLOGY WITH 
INHOMOGENEOUS DARK ENERGY 
EQUATION OF STATE 

In the present section we make brief review of FRW 
cosmology with explicit dark energy equation of state 
(power law). The modification of EOS by Hubble de- 
pendent term (constrained by energy conservation law) 
is done and its role to FRW cosmology evolution is in- 
vestigated. The starting FRW universe metric is: 



ds 2 = -dt 2 + a(t) 2 J2(d^) 2 ■ 



(1) 



In the FRW universe, the energy conservation law can be 
expressed as 



= p + 3i? (p + p) 



(2) 



Here p is energy density, p is pressure. The Hubble rate 
H is defined by H = d/a. When p and p satisfy the 
following simple EOS: 

P = wp , (3) 
and if w is a constant, Eq.J5J can be easily integrated: 

p = Po a-^ +w K (4) 
Using the first FRW equation 
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the well-known solution follows 

a = a (t - ti) 1 *^ or a {t 2 - t)^+^ , (6) 
when w ^ — 1, and 



a e" 



(7) 



when w = — 1. In ©, t\ and t 2 are constants of the in- 
tegration. Eq.Q expresses the deSitter universe. In JBJ, 
since the exponent 2/3(w + 1) is not integer in general, 
we find t > t\ or t < t% so that a should be real number. 
If the exponent 2/3(w + 1) is positive, the first solution 
in Ijfjj expresses the expanding universe but the second 
one expresses the shrinking universe. If the exponent 
2/3(w + 1) is negative, the first solution in {BJ expresses 
the shrinking universe but the second one expresses the 
expanding universe. In the following, we only consider 
the case that the universe is expanding. Then for the 
second solution, however, there appears a singularity in a 
finite time at t = t 2 , which is called the Big Rip singular- 
ity ( for discussion of phantom cosmology near Big Rip 
and related questions, seep], @ and references therein) 
when 



w < — 1 



(8) 



In general, the singularities may behave in differ- 
ent ways. One may classify the future singularities as 
following 0: 

• Type I ("Big Rip") : For t -> t s , a ->• oo, p ->• oo 
and \p\ — ► oo 

• Type II ("sudden") : For t — > t s , a — > a s , p — > p s 
or and \p\ — > oo 

• Type III : For t — > t s , a — > a s , p — > oo and \p\ — * oo 

• Type IV : For t — > t s , a — * a s , p — > 0, \p\ — ► and 
higher derivatives of H diverge. This also includes 
the case when p (p) or both of them tend to some 
finite values while higher derivatives of H diverge. 

Here t s , a s and p s are constants with a s ^ 0. The type 
I may correspond to the Big Rip singularity 1], which 
emerges when w < —1 in J3J). The type II corresponds 
to the sudden future singularity Q at which a and p 
are finite but p diverges. The type III appears for the 
model with p = —p — Ap a 0], which is different from 
the sudden future singularity in the sense that p diverges. 
This type of singularity has been discovered in the model 
of Ref. 3] where the corresponding Lagrangian model of 
a scalar field with potential has been constructed. 
One may start from the dark energy EOS as 



P 



f(p), 



(9) 



^ = P 



(5) 



where /(/?) can be an arbitrary function in general. The 
function /(p) oc p a with a constant a was proposed in 
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Ref. and was investigated in detail in Ref.0]. Using 
iPJ for such choice, the scale factor is given by 



a = ao exp 



dp_\ 



3 J f( P ) J ' 

Using JSJ the cosmological time may be found 

dp 



t = 



K-V3pf(p) ' 



In case 



f(p) = A P a 
by using Eq. l|lU|l . it follows 



a = clq exp 



3(1 - a) A 



(10) 



(11) 



(12) 



(13) 



When a > 1, the scale factor remains finite even if p goes 
to infinity. When a < 1, a — *■ oc (a — » 0) as p — *■ oo for 
A > (A < 0). Since the pressure is now given by 



P : 



P-Ap a , 



(14) 



p always diverges when p becomes infinite. If a > 1, the 
EOS parameter w — p/p also goes to infinity, that is, 
w — > +oo (-co) for A < (A > 0). When a < 1, we 
have to — > -1 + (-1 - 0) for A < {A > 0) as p -> oo. 
By using Eq.lJTIJ for (fT2*|l . one findsQ 



-ct + l/2 



V3kA 1 - 2a 



for a/^, 



and 



In 



t = t 



\/3kA 



for 



(15) 



(16) 



Therefore if a < 1/2, p diverges in an infinite future or 
past. On the other hand, if a > 1/2, the divergence of 
p corresponds to a finite future or past. In case of finite 
future, the singularity could be regarded as a Big Rip or 
type I singularity. 

For the choice (|12J) . the following cases were discussed 
1: 

• In case a = 1/2 or a = 0, there does not appear 
any singularity. 

• In case a > 1, Ea. l|15(l tells that when t — > to, the 
energy density behaves as p — > oo and therefore 
\p\ — * oo due to (|14|) . Ea. (|13fl shows that the scale 
factor a is finite even if p — > oo. Therefore a > 1 
case corresponds to type III singularity. 

• a = 1 case corresponds to the case @ if we replace 
— 1 — A with w. Therefore if A > 0, there occurs 
the Big Rip or type I singularity but if A < 0, there 
does not appear future singularity. 



• In case 1/2 < a < 1, when t — > to, all of p, \p\, and 
a diverge if A > then this corresponds to type I 
singularity. 

• In case < a < 1/2, when t — > to, we find p, 
\p\ — > and a — > ao but by combining (|13|l and 
fT5|l . we find 



ln< 



rT72 



(17) 



Since the exponent (a — l)/(a — 1/2) is not always 
an integer, even if a is finite, the higher derivatives 
of H diverge in general. Therefore this case corre- 
sponds to type IV singularity. 

• In case a < 0, when t — > to, we find p — > 0, a — > ao 
but |p| — * oo. Therefore this case corresponds to 
type II singularity. 

Hence, the brief review of FRW cosmology with specific 
homogeneous EOS as well as its late-time behaviour (sin- 
gularities) is given (see Q for more detail). 

At the next step, we will consider the inhomogeneous 
EOS for dark energy, so that the dependence from Hub- 
ble parameter is included in EOS. The motivation for 
such EOS comes from including of time-dependent bulk 
viscosity in ideal fluid EOS [2( or from the modification of 
gravity (see Appendix). Hence, we suggest the following 
EOS 



P=-p + .f(p)+G(H) . 



(18) 



where G{H) is some function. Then the energy conser- 
vation law J2J has the following form: 



= p + W(f(p) + G(H)) 



(19) 



By using the first FRW equation (JSJ) and assuming the 
expanding universe (H > 0), one finds 

p = F{p)^-^J^( K f{p) + G(n^p~jt)) . (20) 
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G(H) = -f (3H 2 /k 2 ) + —H 



(21) 



Hence, one can express G(H) in terms of / as above. 

As a first example, let assume that EOS© could be 
modified as 



p = wop + w\H 



Using J5J, it follows 



K 2 Wl 



P= \w + 



Therefore w is effectively shifted as 



W — > W e g = W + 



K 2 Wi 



(22) 



(23) 



(24) 
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Then even if wq < — 1, as long as w c s > — 1, there does 
not occur the Big Rip singularity. From another side one 
can start with quintessence value of wq, the inhomoge- 
neous EOS l)23|) with sufficiently negative W\ brings the 
cosmology to phantom era. 

As a second example, we assume f(p) l(T^|) is modified 

as 

f(p) = Ap a - f(p) + G(H) = -Ap a - BH 2fi . (25) 

By using the first FRW equation (J5J, we find f(p) is 
modified as 

MP) = f{ P ) + G(H) = -Ap a - B'pP , 

B '= B (j) ■ ( 26 ) 

If (3 > a, when p is large, the second term in l|26l) becomes 
dominant: 

f eS (p) - B'pP . (27) 

On the other hand, if /3 < a, the second term becomes 
dominant and we obtain l|27|) again when p — ► 0. In 
case of (tl"2l without G(H), when 1/2 < a < 1, there 
is the type I singularity where p goes to infinity in a 
finite time. When G(H) is given by 1)25(1 . if /3 > a, the 
second term in (|26|l becomes dominant and therefore if 
(3 > 1, instead of type I singularity there occurs type 
III singularity. In case of (|12fl with a > 1, the type III 
singularity appears before G(H) is included. Even if we 
include G(H) with /? > a > 1, we obtain the type III 
singularity again and the structure of the singularity is 
not changed qualitatively. For IL'l'l) without G(H), when 
< a < 1/2 or a < 0, there appears the type IV or type 
II singularity where p tends to zero. Since the second 
term becomes dominant if {3 < a, if (3 < 0, the type IV 
singularity for < a < 1/2 case becomes the type II 
singularity but the type II singularity for a < is not 
qualitatively changed. 

In accordance with the previous cases, one finds 

• In case a > 1, for most values of (3, there occurs 
type III singularity. In addition to the type III 
singularity, when < (3 < 1/2, there occurs type 
IV singularity and when (3 < 0, there occurs type 
II singularity. 

• a = 1 case, if (3 > 1, the singularity becomes type 
III. [3=1 case corresponds to (12211 . If (3 < 1 and 
A > 0, there occurs the Big Rip or type I singular- 
ity. In addition to the type I singularity, we have 
type IV singularity when < (3 < 1/2 and type II 
when (3 < 1. 

• In case 1/2 < a < 1, one sees singularity of type III 
for (3 > 1, type I for 1/2 < (3 < 1 (even for (3=1/2) 
or (3 = 1 and B' > (B > 0) case. In addition to 
type I, type IV case occurs for < (3 < 1/2, and 
type II for /3 < 0. 



• In case a = 1/2, we have singularity of type III for 
P > 1, type I for 1/2 < /3 < 1 or /3 = 1 and B' > 
(B > 0), type IV for < /3 < 1/2, and type II for 
(3 < 0. When /3 = 1/2 or (3 = 0, there does not 
appear any singularity. 

• In case < a < 1/2, we find type IV for < (3 < 
1/2, and type II for (3 < 0. In addition to type IV 
singularity, there occurs singularity of type III for 
(3>1, type I for 1/2 < (3 < 1 or (3 = 1 and B' > 
(B > 0) case. 

• In case a < 0, there will always occur type II 
singularity. In addition to type II singularity, we 
have a singularity of type III for (3 > 1, type I for 
1/2 < (3 < 1 or f3 = 1 and £?' > (£? > 0) case. 

Thus, we demonstrated how the modification of EOS 
by Hubble dependent, inhomogeneous term changes the 
structure of singularity in late-time dark energy universe. 

We now consider general case and assume F(p) in H2U|) 
behaves as 

F(p) ~ F oP a , (28) 

with constant Fq and a in a proper limit (e.g. for large p 
or small p). Then when a/1, Ea. l2U|) can be integrated 
as 

F (t - t c ) ~ ^ , (29) 
1 — a 

that is, 

p~ ((l-a)F (*-t c ))^ . (30) 

Here t c is a constant of the integration. When a = 1, the 
energy becomes 

p = Po e Fot , (31) 

with a constant of integration po. By using the first FRW 
equation |[5j), the scale factor may be found 

fl = aoe ±I^f?^(( 1 -^(*-^)^ T , (32) 
when a / 1 and 

a = a e F o v 3 ) (33) 

when a = 1. 

In Q , there has been given an explicit example of the 
EOS where crossing of w = — 1 phantom divide occurs: 

a(t)=a ( J ±-y . (34) 

Here n is a positive constant and < t < t s . The scale 
factor diverges in a finite time (t — > t s ) as in the Big Rip. 
Therefore t s corresponds to the life time of the universe. 
When ( < f S) a(t) evolves as t n , which means that the 
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effective EOS is given by w = — 1 + 2/(3n) > -1. On the 
other hand, when t <~ t g , it appears w = — 1 — 2/(3n) < 
— 1. The solution 134(1 has been obtained with 



4n 

IT 



K 2 p 



(35) 



Therefore the EOS needs to be double- valued in order for 
the transition to occur between the region w < — 1 and 
the region w > — 1. Then in general, there could not be 
one-to-one correspondence between p and p in the above 
EOS. In such a case, instead of l(T%)l . we may suggest the 
implicit, inhomogeneous equation of the state 



F(p,p,H)=0. 
The following example may be of interest: 



(36) 



(37) 



Here Co and Hq are positive constants. Combining l|37|) 
with the energy conservation law 1(19(1 and the first FRW 
equation J5J), one can delete p and p as 



H 2 



!) 



Cnff 4 1 - 



Ho 
H 



which can be integrated as 



H 



1G 



9C 2 H Q (t-t^)(t + -t) 



Here 



t± = t ± 



(38) 



(39) 



(40) 



and to is a constant of the integration. Hence 



3 3 c^ 2 k2 (t-t-y(t + -ty 



(41) 



In l|39|l . since t_ < < t+, as l° n g as t_ <<<<+, 
the Hubble rate H is positive. The Hubble rate H has 
a minimum = Hq when t = to = (t_ + 1 + ) /2 and 
diverges when t — > t±. Then we may regard t — > t_ as a 
Big Bang singularity and t i + as a Big Rip one. As 
clear from l|41(l . the parameter w = p/p is larger than —1 
when t— < t < to and smaller than — 1 when to < t < <+ . 
Therefore there occurs crossing of phantom divide w = 
— 1 when t = to thanks to the effect of inhomogeneous 
term in EOS. 

One more example may be of interest: 



(P + P) + 



16 

kH 



(ho - H) In 



hp-H 
hi 



= . 



(42) 



Here to, ho, hi are constants and ho > hi > 0. A solution 
is given by 



H = hi- hie~ t2/t o , p=^-(hi- hie- t2/t ° 

K \ 

p=-^(hi-hie-^o) 



2\ 2 Ahit _ t 2 , t 



Hence, 



H 



Ihxt 



-t 2 /t 2 o 



(43) 



(44) 



Using the energy conservation law 119|) and the first FRW 
equation JSJ, the second FRW equation may be found: 



2 • 

-— H = p + p 



(45) 



As in 144|) . H is negative when t < and positive when 
t > 0. Eg. (|45[> tells that the effective parameter w — p/p 
of the equation of the state is w > —1 when t < and 
w < — 1 when t > 0. As we find the Hubble rate H 
goes to a constant ho, H — * ho, in the limit of t — > ±oo, 
the universe asymptotically approaches to deSitter phase. 
Therefore there does not appear Big Rip nor Big Bang 
singularity. 

Hence, we presented several examples of inhomoge- 
neous EOS for ideal fluid and demonstrated how the fi- 
nal state of the universe filled with such fluid changes if 
compare with homogeneous case. The ideal fluid with 
implicit EOS may be used to construct the cosmologies 
which cross the phantom divide. 

The interesting remark is in order (see also Appendix) . 
In principle, the more general EOS may contain the 
derivatives of H, like H, H, ... Then more general EOS 
than H36fl has the following form: 



F(p,p,H,H,H,--- ) =0 
Trivial example is that 



p = wp 



2 H 3 ( 1 + U V 



(46) 



(47) 



By using the first JSJ or second (|45(l FRW equations, we 
find 
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P=—H\ p=-—H-—H< 



(48) 



Therefore Eq. (|47(l becomes an identity, which means that 
any cosmology can be a solution if EOS 147(1 is assumed. 
Another, non-trivial example is 



2 • 



p = wp-G t?H + G\H 

K 



>2 



(49) 



Here it is supposed Gq(1 + w) > 0. If Gi(l + w) > 0, 
there appears a solution which describes an oscillating 
universe, 



H = ho cos tot , a = ape 



— ^ sin cut 



(50) 
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Here 



ho = Kt 



Go 



3(1 + w) 



3(1 + w) 

u= \-g^-- (51) 



In case C?i(l + w) < 0, another cosmological solution 
appears 

H = h cosh ut , a = a e^ sinh Qt . (52) 
Here ho is defined by (|51|l again and Co is defined by 



oj = \ -- 



3(1 + w) 



G\K 2 



(53) 



where 5 is a constant. Then combining Eq. I|58[) with 
Eqs. (|56[) . the explicit solution follows 



2 /l 



1 



3 \ t ts-t; 

4/1 1 \ 1 



1 



1 



- - 

3k 2 \t t.-t J t s -t 



where 



5n 2 ' 



(60) 
(61) 
(62) 



(63) 



One can go further and present many more examples of 
inhomogeneous EOS cosmology. 



III. FRW COSMOLOGY WITH 
INHOMOGENEOUS INTERACTING FLUIDS 

In the present section, we study FRW universe filled 
with two interacting fluids. Note that there is some inter- 
est to study the cosmology with homogeneous interacting 
fluids 0, • The inhomogeneous terms for such cosmol- 
ogy may be again motivated by (bulk) viscosity account 

Let us consider a system with two fluids, which satisfy 
the following EOS: 

Pi, 2 = —pi,2 - h,2 (pi, 2) - Gi, 2 {H) . (54) 

For simplicity, the only case is considered that 

P± = w±p± - G± (H) . (55) 

In the above equation and in the following, the indexes 
± instead of 1,2, as pi_2 — P± are used. In a spatially 
flat FRW universe with a scale factor a, the cosmological 
equations are given by 



p± + 3H{p± +p±) = tQ , 

„2 

' P-+P- 

.2 



H = -y(p++P4 



K 



H Z = -{ P+ +P-) 



(56) 
(57) 

(58) 



Not all of the above equations are independent, for ex- 
ample, Eqs.© and © lead to l|57|l. From Eqs.©, 
(|57|l . and the equation for p+ and p+ of H56|) . one obtains 
the equation for p_ and p- of (|56|) .In 4] , the following 
case has been considered: 



In H6Q(I . it is assumed < t < t s . The Hubble rate 
H diverges in a finite time (t — > t s ) as in the Big Rip 
singularity. Therefore t s corresponds to the life time of 
the universe. When t -C t s , H behaves as 2/3t, which 
means that the effective EOS is given by w c s ~ > — 1. 
On the other hand, when t ~ t Sl it appears w c s = —2 < 
— 1. Therefore the crossing of phantom divide w c ff = — 1 
occurs. 

From JSS]-[SB1), we obtain 



± 



1 



W-y. — W 
1 + 



G+(H)+G-{H) 



3 l . w++w- ^ h2 



H 



1 



W+ — W- 



(G' + (H) + G'_{H)) H 



K 

+3H I 1 
1 



w + — w 
1 + 



h 



G+(H) + G-(H) 
2 



3 / . w+ + w_ ^ H - 2 



tH 



9(w + - w-) 3 
4V H 

-\h{G + {H)-G.{H)) 



(64) 



(65) 



First, the case is considered that the Hubble rate H sat- 
isfies the following equation: 



H = S(H) , 



(66) 



G±(H)=0, Q = 5H 2 , w+=0, w_ = -2(59) 



where S(H) is a proper function of H. Hence, Q can be 
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presented as a function of H as 

Q = Q(H) 

' {G' + (H) + G'_(H))S(H) 

H (G + (H) + G-(H)) 



1 

w+ — w. 

W+ + W- 



+3 1 



12 



K 2 — W_) 



1 + 1 HS(H) 



9 



K 2 (w+ — W-) 



X< 1 



K 2 (w+ — W-) 



S'(H)S(H) 



--H{G+(H)-G-(H)) 



(67) 



If Q is given by $7$ for proper G P (H) and S(H), the 
solution of Eqs. (|55l - 1581) can be obtained by solving 
Ea. (|66JI with respect to H. Then from Ij64(l . one finds 
the behavior of p±. As an example, if we consider S(H) 
given by 



S(H) = -—(H-h ) , 
hi 



(68) 



the solution of (|67|l is given by 



the solution of Eas. (jE51 - IS5|l is given by 



H = ho + hi sin (tot + a) 
3 

( h.n -A- hi sin (uit. -4- rv^ 

2/i 2 



p± = — (/io + /iisin(a;t + a)) 



T^5— cos (wi + a) 



(75) 



Here h\ and a are constants of the integration. This 
demonstrates how the inhomogeneous term modifies late- 
time cosmology. 

Choosing G±(H) and Q, one may realize a rather gen- 
eral cosmology. As was shown, if we introduce two fluids, 
even without assuming the non-linear EOS as in 13611 . 
the model crossing w = — 1 effectively can be realized. 
In fact, from l|75[l one has 



H = h\Lo cos {ijjt + a) , 



(76) 



which changes its sign depending on time. When H > 0, 
effectively w < — 1, and when H < 0, w > — 1. Note 
that, as a special case in (|73|) . we may choose, 



G(H) 



3k 2 w 



i _ hi 

H 



S" 2 ■ to 



which gives Q — 0. As Q = 0, from l|56|l . there is no direct 
interaction between two fluids. As is clear from H75JI . 
however, there is an oscillation in the energy densities, 
which may indicate that there is a transfer of the energy 
between the fluids. Hence, the G(H) term might generate 
indirect transfer between two fluids. 



H = h 



t-to' 



(69) 



Here to is a constant of the integration. In the solution 
as H behaves as H ~ j^- when t — t ~ 0, the 

effective w c g is given by w = — 1 + g^-. On the other 
hand, as H becomes a constant ho when t is large, we 
obtain the effective w Q s = 1. 

Next the simpler case is considered: 

w±=-l±w, G±(H)=±G(H). (70) 

Then l|64|) and l|65|l have the following forms: 



P± 



(71) 



2k 2 k 2 w 
1 •• 9u; 

Q = -^H-—H 3 -3HG{H). (72) 

Thus, for example, for an arbitrary G(H), if Q is given 
by a function of H as 



9w 



Q = ^-{H-h )- ^H 3 - 3HG(H) , (73) 



that is, 



H = LU 2 {h -H) , 



(74) 



IV. CROSSING THE PHANTOM BARRIER 

WITH INHOMOGENEOUS EOS AND 
THERMODYNAMIC AL CONSIDERATIONS 

Let us start from the EOS ©. Assuming that w 
crosses —1, which corresponds to f(p) = 0, in order that 
the integrations in (|10fl and (|llf> are finite, f(p) should 
behave as 



hip-poY , o< s <i 



(78) 



Here /(po) = 0. Since < s < 1, f(p) could be multi- 
valued at p = po, in general. Near p = po, Eq. l|ll|) gives, 



- Po) 



l-s 



kV3po/o(1 - s) 
Here t = to when p — po- Since 

# = y/(p), 
from the second FRW Ea. (|45|l . one finds 

ff 2^ /o {-IT 



(79) 



(80) 



1 



k\/3/9 /o(1 - s) 



(81) 
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Hence, when s/(l — s) is positive odd integer, the sign of 
H changes at t = t^, which shows the crossing w = — 1. 

In recent paper based on consideration of mixture 
of two fluids: effective quintessence and effective phan- 
tom, the following, quite interesting EOS has been sug- 
gested: 



Ap m + Bp m = [Cp m + Dp m ) c 



(82) 



Here A, B, C, D, and a are constants and m is an integer. 
This EOS can be regarded as a special case of By 
writing p as 

p = Q( P )p , (83) 

one obtains 

^m(a-i) _ p (Q m ) = (A + BQ m ) (C + DQ m y a . (84) 
Since 

-a-l 



F'{Q m ) = (C + DQ r ' 

x (BC - aAD + (1 - a)BDQ r 

it follows F' (Q m ) = when 
Q m - 



D B 

1 — a 



By properly choosing the parameters, we assume 

1 . 



D B 



I - a 



(85) 



(86) 



(87) 



When Q 



1. 



Here 
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F(Q m ) ~q Q + q 2 (Q + IT 



F(l) = {A-B)(C + D)~ a , 

1 cPF_ 

2 dQ* Q= _ x 



-a{a-l){C - D)- a ~ z D\A- B)m z . (89) 



In l|89(l . it is supposed m is an odd integer. Solving 184(1 
with l|88|l with respect to Q, one arrives at 



Q 



-l ± 



m(a — l)pg 



ma— m — 1 



(p-po) 



'12 



1/2 



Here po is defined by 



9o = Pq 

Using JS3J), the function Q is 



(a-l) 



m(a - 



ma— m — 1 



(P - Po) 



<72 



1/2 



(90) 



(91) 



(92) 



Comparing (§2J with (JTHJ, we find that the EOS (52jl 

surely corresponds to s = 1/2 case in (I78f) . 

For the EOS (|82jl . there are interesting, exactly solv- 
able cases. We now consider such a case and see that 
there are really the cases of EOS crossing barrier w = — 1. 
The energy conservation law J5J may be rewritten as fol- 
lows: 



p = -p-V%, V^Voa* 
dv 



(93) 



Here Vq is a constant with the dimension of the volume. 
Use of Eq. (JHSl gives 



= V^ + (1 + Q(p))p 



(94) 



Using l|84|l . we further rewrite (|94|) as an equation with 
respect to Q: 

(BC - aAD + BP (I — a) Q m ) Q 171 ' 1 dQ 
" ~ (l-a)(C + DQ m ){A + BQ m ) dV 

+1 + Q. (95) 

Assuming Eq. Q87JI , the above Eq. Q95J1 takes a simple form: 







BDil + Q^Q^ _ v dQ +l + Q . (96) 



(C + L>Q m ) (A + BQ m ) dV 



Especially in the simplest case m = 1, one can easily 
solve pjj) 



Q 



GO' 



D - B 



(*) 



/3 



Here Vi is a constant of the integration and 

I 



/3 = 



AD-BC (1-a) |-1 



(97) 



(98) 



In the above equation, Ea. H87fl is used. Hence, when 
(U/Vi)' 3 -> 0, it follows w = p/p = Q -► -C/L>. On 
the other hand, when (V/Vi) — > oo, one arrives at W = 
Q — > —A/B. Hence, the value of w changes depending 
on the size of the universe. Especially when 



V (C -D 



Vi \A-B 



1//3 



(99) 



there occurs the crossing of phantom divide w = Q = — 1 
(compare with 01). 

As the inhomogeneous generalization of the EOS l|82jl. 
we may consider 



A I —H' 



Bp m = (Cp m + Df 



or 



Ap m + Bp m = [C[ —H< 



Dp r ' 



(100) 



(101) 
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or, more general EOS 



As an example, for the solution (|75|l 



(A - A')p m + A 1 ( ^H' 



= [{c- c) P m + c 4 ff 



+ Bp m 
Dp m \ . (102) 



By using the first FRW equation ©, it folows that the 
EOS (Pb|) . fTUT|) . and (JT^J) are equivalent to Es- 
pecially if to = 1 and (|87|l could be satisfied, one obtains 
the solution l(9"7|) . 

Hence, using the first and second FRW Eqs.JSJ) and 
(gnj, the EOS (|H2J) with m = 1 can be rewritten as 



dt 2 



a 2[ - B > 



3k 2 (A-B) /3k 2 (C-D) 



AB 2 



AD 2 



3 

Xtt 2 



{^-^-%)}^{a^-%))\ .CU).-!! 



When (|87|l is satisfied, this second order differential Eq. 
looks as 



AB 



d 2 X 

~d£ r \Zk 2 (A-B) 



a-l 



d 2 Xl 
dt 2 



(104) 



which also admits, besides the solution crossing w = — 1 
(|97jl . a flat universe solution 



a = clq , (ao : constant) , 
and deSitter universe solution 



a = aoe 



(105) 



(106) 



As next generalization of 1)82(1 . one may consider the 
following EOS: 



Ap + Bp - 



A-B 



-H 2 



Cp + Dp-^-^-H 2 



(107) 



Here a is assumed to be a function of H. Then by using 
the first and second FRW Eqs.JSJ and gU, the EOS 
(|107fl can be rewritten as 



2D h 



a (H) 



which gives 



(108) 



(109) 



= — 
hi 



H 



dH 



(110) 



1 - 



H-hr 



Comparing l|109|) with H110|) . in case that 



2D 



, a{H) = l 



2 ki- 



ln 1 



, (HI) 



the solution J^HJ) follows from the EOS ifTHTjl . 

As another generalization of 1(82(1 , we may consider the 
following EOS: 



Ap m + Bp m = G(H) (Cp m + Dp" 



(112) 



Here G(H) is a function of the Hubble rate. For simplic- 
ity, the following case is considered 



to = 1 



(113) 



Then, writing p as (l8"3")l and using Q, the energy looks 
like 



p = (A + BQ)^+^t(C + DQY 



(114) 



which corresponds to l(84() . Assuming Eg. ((87(1 . by using 
((83(1 . instead of 1(96(1 . one gets 







(1 -a)BD 



(l-a--r)(C + DQ)(A + BQ) dV 
which can be solved as 

J 



^S + l, (H5) 



C - A 



(*) 



Here V\ is again a constant of the integration and 



(116) 



/3: 



(1 - a)BD 



{l-a-i)AD-BC ' 



(117) 



Then as in QJJ, when (V/V^ 



0, we have w 



pi p = Q —> —C/D and when (V/Vi) — > oo, we have 
w = Q — > —A/B. The power of V, however, is changed 
in Eg. 1(116(1 if compare with Eg. ((97(1 . Thus, we presented 
number of FRW cosmologies (including oscillating uni- 
verses) filled by cosmic fluid with inhomogeneous EOS 
where phantom divide is crossing. Definitely, one can 
suggest more examples or try to fit the astrophysical data 
with more precise model of above sort. 

In [l^, t ne thermodynamical models of the dark en- 
ergy have been constructed. Especially it has been shown 
that, for the fluid with constant w, the free energy 
F(T, V) is generally given by 

F(T,V) =TP({T/T ) 1 ' W {V/V Q )) . (118) 
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Here T is the temperature and V is the volume of the 
universe. For the dimensional reasons, the positive pa- 
rameters To and Vq are introduced. 

The interesting question is what happens with the en- 
tropy when the value of w crosses —1. As a model, the 
case that w — Q depends on V as in (|97|) may be consid- 
ered: 



w = w(V) = 



wo + Wl 



(*) 



/3 



1 



(ft)' 



(119) 



When [3 > 0, w — > wq for small universe and «)-t«)i for 
large universe. 

The specific dependence of free energy may be taken 
as below 



F 



f T \f T \ ^vl V 



Vn 



(120) 



Here 7 is a constant. When 7 = 1 and w is a constant, 
the free energy is proportional to the volume. For usual 
matter, due to self-interaction and related effects, 7 is not 
always unity. Then, the pressure p, the energy density p, 
and the entropy iS are given by 



P 



7-1 



s = -— 



dF 
dV 

M (T_ 
V I To 



V \ dT 



fox (t\ 1+ ^ (vy 7 - 1 

wV \T J \V 0/ 
dF 
df 

fa f, , 7\ / T 



To 



1 



w J \T< 



(121) 



In the pressure p, the second term in large { } comes from 
V dependence of w in H119(l , which vanishes for large or 
small universe (V — > 00 or V — > 0). Hence, for small or 
large universe p/p — > w(V) — > wo,i- As seen from the 
expression for S, the sign of the entropy changes at 



w = —7 



(122) 



If 7 = 1 , the sign of the entropy S changes when crossing 
w = — 1 (the entropy becomes negative when w is less 
than —1 as it was observed in 0), but in the case that 



7 < \w \, I tOi| 



(123) 



the entropy does not change its sign. 

We should note that the expressions (|121|l are not 
wcll-dchncd, unless 7 = 0, when w = 0, which corre- 
sponds to dust. One may assume < 7 < wq <ti 1 and 
wi ;$ — 1. Then as clear from l|119|) . w changes from 
Wo ~ for small universe to Wi < — 1 for large universe 
and crosses — 1. Since we always have (7/100 1 < 1 an d 
therefore 1+7/10 > 0, the entropy S (|121fl is always 
positive and does not change its sign as long as /o < 0. 
This explicitly demonstrates very beatiful phenomenon: 
there exist thermodynamical models for dark energy with 
crossing of phantom divide. Despite the preliminary ex- 
pectations, the entropy of such dark energy universe even 
in its phantom phase may be positive! 



V. DISCUSSION 

In summary, the effect of modification of general EOS 
of dark energy ideal fluid by the insertion of inhomoge- 
neous, Hubble parameter dependent term in the late-time 
universe is considered. Several explicit examples of such 
term which is motivated by time-dependent bulk viscos- 
ity or deviations from general relativity are considered. 
The corresponding late-time FRW cosmology (mainly, in 
its phantom epoch) is described. It is demonstrated how 
the structure of future singularity is changed thanks to 
generalization of dark energy EOS. The number of FRW 
cosmologies admitting the crossing of phantom barrier 
are presented. The inhomogeneous term in EOS helps to 
realize such a transition in a more natural way. 

It is interesting that in the case when universe is filled 
with two interacting fluids (for instance, dark energy and 
dark matter) the Hubble parameter dependent term may 
effectively absorb the coupling between the fluids. Again, 
in case of two dark fluids the phantom epoch with possi- 
bility of crossing of w — — 1 barrier occurs is constructed. 
It is also very interesting that there exists thermody- 
namical dark energy model where despite the preliminary 
expectations |17| the entropy in phantom epoch may be 
positive. This is caused by crossing of phantom barrier. 

As it was demonstrated making the dark energy EOS 
more general, this extra freedom in inhomogeneous term 
brings a number of new possibilities to construct the late- 
time universe. One can go even further , assuming that 
inhomogeneous terms in EOS are not restricted by energy 
conservation law (as it is often the case in braneworld 
approach). Nevertheless, only more precise astrophysical 
data will help to understand which of number of EOS of 
the universe under consideration (in other words, dark 
energy models) is realistic. 
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In the situation where the derivative of H can be ne- 
glected as H < H 2 or H < H 3 , we find 



APPENDIX A: INHOMOGENEOUS TERMS 
FROM MODIFIED GRAVITY 

Let us consider the possibility to obtain the inhomoge- 
neous EOS from the modified gravity. As an illustrative 
example, the following action is considered: 



1 



S = / (Pxy/=g ( ——^ + C 



2k 2 



f(R) 



(Al) 



Here f(R) can be an arbitrary function of the scalar cur- 
vature R and £ ma tter is the Lagrangian for the matter. 
In the FRW universe, the gravitational equations are: 







H+p— f[R = 6H + 12H- 



dt 



6 [H + H 2 - H 



xf (R = 6H+12H : 



(A2) 



2H + 3H 2 ) +p + f (R = QH+ 12 H 



2 -H - 3H Z 



dt 2 

xf [R= 6H+12H 2 



dt 



(A3) 



Here p and p are the energy density and the pressure 
coming from £ ma tter- They may satisfy the equation of 
state like p — wp. One may now define the effective 
energy density p and p by 



p-f[R= 6H+12H' 



6 [ H + H — H 



dt 



xf [R = 6H+12H : 



(A4) 



p + f ( R = 6H + 12H- 
d 2 



2 —H — 3H 



dt 2 



dt 



x f (r = 6i? + 12H 2 } . 
Thus, it follows 

p = wp+ (1 + w)f (R = 6H + 12H 2 



(A5) 



2 I (-1 - 3w) H - 3 (1 + w) H z + 
(2 + 2>w) f (R = 6H + 12H 2 ^ .(A6) 



p ~ wp + G(H) , 
G{H) = (l + w)f(R= 12H 2 ) 

-3(l + w)H 2 f (R = 12H 2 ) 



(A7) 



Typically H has a form like H ~ ho/ (t — tx) or H ~ 
ho/ (*2 — t), with ho = 2/3 (w + 1), corresponding to J()J). 
Hence, the condition H <C H 2 or H <C H 3 requires ho ^S> 
1, which shows w ~ — 1 as in the modern universe. This 
supports our observation that inhomogeneous terms may 
be the effective ones which are predicted due to currently 
modified gravity theory. 

The modification of the EOS by G(H) terms might 
come also from the braneworld scenario. Indeed, the sin- 
gle brane model is described by the following simple ac- 
tion 



M 2 

S = ^ 



d*xdyy/-glftR& 



J d*V=$ (M^R + Anattor) 



(A8) 



Here Mpj = l/8wG, y is the coordinate of the extra di- 
mension, and £ ma ttcr is the Lagrangian density of the 
matters on the brane. The five-dimensional quantities 
arc denoted by suffix "(5)". Inref.[3 it has been shown 
that the FRW equation for 4d brane universe could be 
given by 



(A9) 



Here p is the matter energy density coming from £ ma tter- 
More general case is considered in ref.^| where the FRW 
equation is modified as 



3 R 2 tt a 

K 2 \ 



(A10) 



Here a is a constant. One may assume that the matter 
energy density p satisfies the energy conservation as in 
©. Then from JXSJ, wc find 

-^('-ipr) *-' + '• < A11 » 

By comparing l)A10jl with the first FRW equation (0 
and ijAllfl with the second FRW equation l|45|). one may 
define the effective energy density p and pressure p as 



p = p- 



3H a 
^P 7 



3H C 



K 2 r 2 c - a 



aH a - 2 H 
n 2 r 2 c - a 



(A12) 



They satisfy the first JSJ and second (|43|) FRW equations: 



K 



;H = p + p 



(A13) 
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If it is also assumed the matter energy density p and 
the matter pressure p satisfy the EOS like p = wp, the 
effective EOS for p and p is given by 

p = wp-(l + w) . (A14) 



This shows that brane-world scenario may also suggest 
various forms of inhomogencous modification for effective 
EOS of matter on the brane. 



Especially if one can neglect H, it follows 

p~wp-(l + w) 2 2 _ a . (A15) 

K T c 
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